Exchange and correlation effects in the relaxation of hot electrons in noble metals 
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We report extensive first-principles calculations of the inelastic lifetime of low-energy electrons in 
the noble metals Cu, Ag, and Au. The quasiparticle self-energy is computed with full inclusion of 
exchange and correlation (xc) effects, in the framework of the GWT approximation of many-body 
theory. Although exchange and correlation may considerably reduce both the screening and the 
bare interaction of hot electrons with the Fermi gas, these corrections have opposite signs. Our 
results indicate that the overall effect of short-range xc is small and GWT linewidths are close to 
their xc-free G°W° counterparts, as occurs in the case of a free-electron gas. 
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I. INTRODUCTION 

Relaxation lifetimes of excited electrons with ener- 
gies larger than 1 eV are mainly dominated by electron- 
electron (e-e) inelastic interactions of the excited elec- 
tron with the electrons in the valence bands of the solid. 
The inelastic lifetime of these so-called hot electrons has 
been investigated for many years on the basis of the free- 
electron gas (FEG) or jellium description of the solidji* 2 * 3 - 
in which a homogeneous assembly of interacting electrons 
is assumed to be immersed in a uniform positive back- 
ground. Nonetheless, time-resolved two-photon photoe- 
mission (TR-2PPE) experimenta 4 i 5 i 6 i 7 i 8 iV°i 1:L and ballis- 
tic electron-emission spectroscopy (BEES)^ 2 " have shown 
that band-structure effects play a key role in the decay 
mechanism. 

The first theoretical investigations of the hot-electron 
dynamics which take into account explicitly the band 
structure of the solid were performed only a few years 
a g 0> i3 i i4 g mce then, first-principles calculations of the e- 
e scattering have been reported for simple^ noble 
and transition metals>i2i22i2i Nevertheless, all existing 
calculations have been performed within the G°W° ap- 
proximation of many-body theory^ with no inclusion of 
exchange and correlation (xc) effects. 

In this work, we go beyond the G°W° approximation 
by including xc effects both in the description of the dy- 
namical screening of the many-electron system [we go 
beyond the random-phase approximation (RPA) in the 
evaluation of the screened interaction W] and in the ex- 
pansion of the electron self-energy in terms of W. This is 
the GWT approximation* 2 ^ 2 ^ which treats on the same 
footing xc effects between pairs of electrons within the 
Fermi sea (screening electrons) and between the hot ex- 
cited electron and the Fermi sea. Mahan and Sernelius 23 
showed that the inclusion of the same vertex function 
in the screened interaction and the numerator of the self- 
energy yields results for the bandwidth of a homogeneous 



electron gas very similar to those obtained in the G a W° 
approximation, due to a large cancellation of vertex cor- 
rections. Nonetheless, due to the well-known differences 
between the dynamical response of the noble metals and 
the homogeneous electron gas, the impact of exchange 
and correlation on the relaxation of hot electrons in these 
materials required further investigation. 

Our first-principles calculations start by solving self- 
consistently the Kohn-Sham equations of density func- 
tional theory (DFT)f2£*2& within the local density ap- 
proximation (LDA) and with the use of a plane-wave ex- 
pansion of the single-particle Kohn-Sham Bloch states. 
The electron-ion interaction is described by a nonlocal, 
norm-conserving ionic pseudopotential, 27 by keeping all d 
electrons as valence electrons. The single-particle Bloch 
states are then used to compute the GWT electron self- 
energy and hot-electron inelastic lifetimes. For compar- 
ison, we also compute G°W°, G°W, and GW°T decay 
rates, with no inclusion of xc effects, with inclusion of xc 
effects beyond the RPA in the screened interaction W, 
and with inclusion of xc effects beyond the G°W° in the 
expansion of the electron self-energy in terms of the RPA 
screened interaction W°, respectively. Our results indi- 
cate that the overall effect of short-range xc is small and 
GWT linewidths are close to their G°W° counterparts, 
as occurs in the case of a FEG. 

In order to establish the role that occupied d states 
play in the relaxation of hot electrons in the noble met- 
als, we also use a pseudopotential with all d electrons 
assigned to the core. We find that a major contribution 
from occupied d states participating in the screening of 
e-e interactions yields lifetimes of low-energy excited hot 
electrons that are larger than in the absence of d states 
by a factor of ~ 2 in Cu and by a factor of ~ 2.5 in Ag 
and Au. 

This paper is organized as follows. Explicit expressions 
for the electron decay rate (inelastic lifetime broadening) 
in a FEG and periodic crystals are derived in Sec. [H] in 
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the GWT approximation of many-body theory. The re- 
sults of numerical calculations of hot-electron lifetimes in 
the noble metals Cu, Ag, and Au are presented in Sec. lIIII 
The summary and conclusions are given in Sec. IIVI Un- 
less stated otherwise, atomic units are used throughout, 



1. e. 



= H = m e = 1. 



II. THEORY 

Let us consider an arbitrary many-electron system of 
density no(r). In the framework of many-body theory, 
the damping rate or reciprocal lifetime of a quasiparticle 
in the single-particle state <fii(r) of energy Si (e, > £_f, £f 
being the Fermi energy) is obtained as the projection of 
the imaginary part of the electron self-energy X(r,r';£j) 
over the quasiparticle-state itself^ 



-2 Jdr Jdr' <p*(r) ImS(r, r'; £< )^(r'). 



(1) 



Within many-body perturbation theory, it is possible 
to obtain the electron self-energy as a series in the bare 
Coulomb interaction v(r, r'), but due to the long range of 
this interaction such a perturbation series contains diver- 
gent contributions. Therefore, the electron self-energy is 
usually rewritten as a series in the frequency-dependent 
screened interaction W(r, r'; lu). To lowest order in the 
screened interaction, the self-energy is obtained by in- 
tegrating the product of the interacting Green function 
G(r, r',£i — ui) and the screened interaction W(r, r'; u>), 
and is therefore called the GW self-energy. If one further 
replaces the interacting Green function by its noninter- 
acting counterpart G°(r,r',£i — u), one finds the G°W 
self-energy and from Eq. the following expression for 
the G°W lifetime broadening: 



dr 



dr' </>*(!•) </>}(!■') 



x ImW(r,r';£i - £/) <^(r') 0/(r), 



(2) 



where the sum is extended over a complete set of single- 
particle states 0/(r) of energy ej (ef < £/ < £i)- The 
screened interaction W(y,y' \co) can be rigorously ex- 
pressed as follows 



W(r,r';u) = v(r,r') 



J dri J dr 2 v(r,r 1 



x x(ri,r2;w)u(r 2 ,r'), 



(3) 



x(j,r;uj) being the time-ordered density-response func- 
tion of the many-electron system, which for the pos- 
itive frequencies (cu > 0) entering Eq. J3J) coincides 
with the retarded density-response function of linear- 
response theory. In the framework of time-dependent 
DFT (TDDFT)fS the exact retarded density-response 
function is obtained by solving the following integral 
equation 
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X(r,r';w) = x°(r,r';w) + / dri / dr 2 x°(r,ri;w) 



x {«(n,r 2 ) + r c [no](ri,r 2 ;w)}x(r 2 ,r';cj), (4) 

where x ( r > r 'jk') denotes the density-response function 
of noninteracting electrons 



X °(r,r';c) = 2^ 



fi fj 



IT] 



x &(r)^(r)^(r')&V)- (5) 

Here, 4>i(r) and £j denote the eigenfunctions and eigenval- 
ues of the Kohn-Sham hamiltonian of DFT, /j are Fermi- 
Dirac occupation factors, rj is a positive infinitesimal, and 
the frequency-dependent xc kernel f xc [no](r, r'; us) is the 
functional derivative of the frequency-dependent xc po- 
tential T4 c [n](r,cj) of TDDFT, to be evaluated at no(r): 

SVn C [n](r,u) 



f xc [n }(r 7 r';uj) = 



8n(r',u) 



(6) 



In the RPA, f xc [no](r,r' ;tu) is set equal to zero and 
Eq. yields the so-called G°W° lifetime broadening^ 
In the adiabatic LDA (ALDA)^ 



r c (r,r';^ = 



dV xc (n) 



dn 



5(t-t>), (7) 



=io (r) 



V xc (n) being the static xc potential of a uniform electron 
gas of density n. 

The xc kernel f xc [no](r, r'; u>), which is absent in the 
RPA, accounts for the presence of an xc hole associ- 
ated to all electrons in the Fermi sea. Hence, one might 
be tempted to conclude that the full G°W approxima- 
tion [with the formally exact screened interaction W of 
Eq. should be a better approximation than its G°W° 
counterpart [with the screened interaction W evaluated 
in the RPA]. However, the xc hole associated to the ex- 
cited hot electron is still absent in the G°W approxima- 
tion. Therefore, if one goes beyond RPA in the descrip- 
tion of W, one should also go beyond the G°W approx- 
imation in the expansion of the electron self-energy in 
powers of W. By including xc effects both beyond RPA 
in the description of W and beyond G°W in the descrip- 
tion of the self-energy^2i24 the so-called GWT approxi- 
mation yields a lifetime broadening that is of the G°W 
form [see Eq. but with the actual screened inter- 
action W(t,t';oj) replaced by a new effective screened 
interaction 



W(r, r'; to) = v(r, r') + / dri / dr 2 {v(r, ri) 

+ r> ](r,r i;w )} x(ri,r 2 ; W )t;(r 2 ,r'), (8) 
which includes all powers in W beyond the G a W approx- 



imation. 



A. Free-electron gas 

In the case of a uniform FEG, there is translational 
invariance in all directions, the single-particle states en- 
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tering Eqs. and (JSJ) are momentum eigenfunctions 
</>k( r ) — exp(ik • r) of energy — k 2 /2, and Eq. (J2J 
is easily found to yield 



-2 



(2^ 



Im W(q, £ k - £k-q), 



(9) 



the integral being subject to the condition that Ef < 
£k- q < £k, and W(q, u>) being the Fourier transform of 
the screened interaction W(r,r';u>) of Eq. |JSJ, which in 
the GWT approximation should be replaced by the effec- 
tive screened interaction W(r, t';uj) of Eq. (JSJ. We note 
that Eq. (0) with W(q, u>) replaced by the Fourier trans- 
form of W(t,t';ui) yields precisely the decay rate that 
one would obtain from the GWT self-energy of Ref. I23I 



B. Periodic crystals 

For periodic crystals, the single-particle states enter- 
ing Eq. 0) are Bloch states 0k,i(r) and ^>k- q ,/(r) with 
energies £k,j and £k- q ,/, i and / denoting band indices. 
Hence, Eq. 0) yields 

J = iE / dt * £ £k,i ; k- q ,/(G)B k , i: k- q ,/(G') 

(10) 



f G,G' 
X Im?G,G'(q,Ek,i - e k-q,/), 



where the integral is extended over the first Brillouin 
Zone (BZ), the vectors G and G' are reciprocal-lattice 
vectors, 

£ M;k - q ,/(G) = |dr0 kiJ (r)e^+ G ) r k - q ,/(r), (11) 

and WG,G'(q, w ) denote the Fourier coefficients of the 
screened interaction W(r, r';u), which are usually ex- 
pressed as follows 



^G,G'(q,w) = e G 1 G ,(q,w)u G -(q) 



(12) 



«G(q) = 47r/|q+ G| 2 being the Fourier transform of the 
bare Coulomb interaction, and e G 1 G /(q, w) being the so- 
called inverse dielectric matrix. 

We remark that for a given hot-electron energy e there 
are in general various possible wave vectors and bands. 
Since r^Ji = t^, S denoting a point-group symmetry 
operation in the crystal, one only needs to consider states 
inside the irreducible wedge of the BZ (IBZ). An energy- 
dependent reciprocal lifetime r _1 (e) can then be defined 
by doing a weighed average over all wave vectors and 
bands lying with the same energy in the IBZ. 



1. G°W° approximation 

In the G°W° approximation, the xc kernel 
f xc [no)(r, r'; u>) entering Eqs. (0) and JSJ is set equal to 



zero. Hence, in this approximation the dielectric matrix 
is 

eG,G'(q>^) =5 G ,G> - "G(q)XG,G'(q> w )i ( 13 ) 

Xq g , (q, u>) being the Fourier coefficients of the nonin- 
teracting density-response function of Eq. 0) : 



/k,n /k+q/, 



£k,n - £k+ q ,n' +UJ + ir] 



V«|e~ l(q+G) - r |0k +q ,„'> 



'k+ q ,n'| e 



: (q+G')-r| 



(14) 



Couplings of the wave vector q + G to wave vectors 
q + G' with G ^ G' appear as a consequence of the 
existence of electron-density variations in the solid. If 
these terms, representing the so-called crystalline local- 
field effects (LFE), are neglected, one can write Eq. i|l(J|) 
as follows 



|-Bk,i;k-q,/(G)|' 



^ u D ^ G |q + 
Im [e G ,G(q, £k.» - £k- q ,/)] 



(15) 



|£G,G(q,£k,i - £k-q,/)| 2 

This expression accounts explicitly for the three main in- 
gredients entering the hot-electron decay process. First 
of all, the coupling of the hot electron with available 
states above the Fermi level is dictated by the matrix 
elements i?k,i ; k- q ,/(G). Secondly, the imaginary part of 
the dielectric matrix £G,G(q, £k,i — £k- q ./) represents a 
measure of the number of states available for the cre- 
ation of e-h pairs with momentum and energy q + G and 
£k,i — £k-q,/, respectively. Thirdly, the dielectric matrix 
in the denominator accounts for the many-body e-e in- 
teractions in the Fermi sea, which dynamically screen the 
interaction with the external hot electron. 



2. G°W approximation 

In the G°W approximation, the Fourier coefficients 
WG,G'(q,w) still take the form of Eq. but with 

the inverse dielectric matrix now given by the following 
general expression: 



£ g!c (<!' w ) = s g,g> + uc(q) xg,g> (q, w), 



(16) 



where xg.c (q, are the Fourier coefficients of the in- 
teracting density-response function of Eq. 101 : 

XG,c (q, w) = x°g,g> (q> w ) + £ ^G.d (q 5 w ) 
x {«G 1 (q) + / G < ;,G 2 N(q,w)}xG 3 ,G'(q ! w) ! (17) 

/ G c G z [no] (q, bS) being the Fourier coefficients of the xc 
kernel of Eq. ©. In the ALDA, 

dV xc (n) 



dn 



(18) 



n=n (r) 
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3. GWT approximation 

In the GWT approximation, the lifetime broadening 
is still of the form of Eqs. (JTDJ-dJ, but with the test- 
charge-test-charge inverse dielectric matrix of Eq. 11 tip 
replaced by the new test-charge-electron inverse dielec- 
tric matrix 

G" 

+ /G°G«[«o](q,w)} XG»M<b")- (19) 

where XG,G'(qj ar e the Fourier coefficients of Eq. 1)170. 
If one sets the xc kernel /q c G / [«o] (q, w) entering Eqs. I|17|) 
and i|19|) equal to zero, one finds the G°W° approx- 
imation. Instead, if one only sets the xc kernel 
fee i n o] (qj w ) entering Eq. l(T7|l equal to zero, one finds 
the GW°T approximation. 

III. RESULTS AND DISCUSSION 

In the calculations presented in this section, all 
the single-particle Bloch states and energies entering 
Eqs. I |10| l. (|llfl . and (|14|l are taken to be the eigenfunc- 
tions and eigenvalues of the LDA Kohn-Sham hamil- 
tonian of DFT. We expand the single-particle Bloch 
states in a plane-wave basis, we invoke the LDA 
with the Perdew-Zunger parametrization^fi of the Quan- 
tum Monte Carlo uniform-gas xc energy of Ceperly 
and Alder—, and we describe the electron-ion interac- 
tion with the use of a nonlocal, norm-conserving ionic 
pseudopotential, 27 

Cu, Ag, and Au are noble metals with entirely filled nd 
bands, n being 3, 4, and 5, respectively. Slightly below 
the Fermi level [at e — £f ~ 2 eV in Cu and Au, and at e — 
ef ~ 4eV in Ag] there are d bands capable of holding 10 
electrons per atom. The one remaining (n + l)s electron 
per atom occupies a free-electron-like band below and 
above the d bands. Hence, a combined description of 
the localized nd 10 and delocalized (n + l)s 1 electrons is 
needed in order to address the actual electronic response 
of these materials. 

The results presented below have been found by ei- 
ther keeping all nd 10 and (n -I- l)s 1 electrons (full cal- 
culation) or keeping only the (n + l)s 1 electrons (s- 
calculation) as valence electrons in the generation of 
the pseudopotential. The full calculation has required a 
kinetic-energy cut-off as large as 75 Ry. 37 Well-converged 
results have been found for all hot-electron energies un- 
der study (0.5 eV < e — £f < 3.5 eV), with the inclusion 
of conduction bands up to a maximum energy of ~ 25 eV 
above the Fermi level. Samplings of the BZ have been 
performed on a 20x20x20 mesh containing 256 points in 
the IBZ, although well-converged results have sometimes 
been obtained on 16x16x16 meshes. Reciprocal-space 
sums have been extended over 15 G vectors, the mag- 
nitude of the maximum momentum transfer q+G being 




FIG. 1: Average lifetime of hot electrons in Cu, as a function 
of the hot-electron energy e — ef with respect to the Fermi 
level. Solid and open circles represent full GWT and G°W° 
calculations, as obtained from Eqs. 11011 - 112(1 with the test- 
charge-electron inverse dielectric matrix of Eq. 119H and with 
the RPA dielectric matrix of Eq. 1131 , respectively. Solid and 
open squares represent the corresponding GWT and G°W° 
s-calculations, respectively, where the 3d-shell is assigned to 
the core in the generation of the pseudopotential. Solid and 
dashed lines represent GWT and G°W° FEG calculations, 
as obtained from Eq. JUJ with the electron-density parame- 
ter r a = 2.67 [r s = (3/ATrno) 1 ' 3 ] corresponding to the av- 
erage density no of 4s 1 electrons in Cu. The Fourier coef- 
ficients /g°g' [ n o] (li w ) an d f xc [no](ci, w) entering the crystal 
and FEG calculations have both been calculated in the ALDA 
with use of the Perdew-Zunger xc potential of a uniform elec- 
tron gas. 

well over the upper limit of ~ 2qp , qF being the Fermi 
momentum. 



A. Copper 

Fig. d shows our full first-principles GIFT calculation 
of the average lifetime t(E) of hot electrons in Cu (solid 
circles), as obtained in the ALDA from Eqs. l(TU j) -lJT2 )l but 
with the inverse dielectric matrix £Q G ,(q, u) replaced 
by that of Eq. i|19|) . For comparison, we also plot in 
the same figure our full first-principles G°W° calculation 
(open circles), which reproduces previous calculations, 13 
first-principles GWT and G°W° s-calculations (solid and 
open squares), and GWT and G°W° FEG calculations 
with the electron density no equal to that of 4s 1 electrons 
in Cu (solid and dashed lines). 

We note from Fig. Q that xc effects yield hot-electron 
lifetimes that are larger than in the absence of exchange 
and correlation by no more than ~ 3% in the whole 
hot-electron energy range under study, as occurs when d 
electrons are assigned to the core (s-calculations) and in 
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FIG. 2: Impact of exchange and correlation on the lifetime of 
hot electrons in Cu, as a function of the hot-electron energy 
e — ef wit respect to the Fermi level. Triangles, diamonds and 
circles represent the ratios t g o w /t g o w o , t gw o v /t g o w o , and 
tgwt /t g q w q in the presence (upper panel) and in the absence 
(lower panel) of 3d 10 occupied states, respectively. Dashed, 
dashed-dotted and solid lines represent the corresponding ra- 
tios for a FEG with r a = 2.67. 



the case of a FEG. Exchange-correlation effects included 
in the GWT scheme have two sources, as discussed in 
Sec. [n] Firstly, there is the reduction of the screening 
due to the presence of an xc hole associated to occupied 
states below the Fermi level, as in the G°W approxi- 
mation. Secondly, there is the xc hole associated to the 
excited hot electron over the Fermi level, as in the GW°T 
approximation. These contributions have opposite signs 
and it is the latter which dominates. 

The impact of xc effects on the lifetime of hot electrons 
in Cu is illustrated in Fig. where G°W, GW°T, and 
GWT calculations are compared to their G°W° coun- 
terparts, with (upper panel) and without (lower panel) 
inclusion of d states. The impact of exchange and corre- 
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FIG. 3: Scaled lifetimes r(s) x (e — Ef) 2 of hot electrons in 
Cu. Open circles, squares, triangles, and diamonds represent 
G°W° calculations obtained from Eq. 11511 (hence, with no 
inclusion of crystalline local- field effects), (i) in the presence 
of 3d states (open circles), (ii) in the absence of 3d states 
(squares), (iii) including 3d states only in the evaluation of the 
imaginary part of £G,G(q, w) (triangles), and (iv) including 3d 
states in the evaluation of the imaginary part of £G,G(q, 
as well as in the evaluation of the denominator |eG,G(q, w)| 
(diamonds). Solid circles represent the full GWT matricial 
calculation of Fig. but now multiplied by (e — £f) 2 - The 
dashed line represents the G°W° calculation for hot electrons 
in a FEG with r s = 2.67. 



lation on either screening or the bare interaction of the 
hot electron with the Fermi sea, which is significant in 
the absence of d states (lower panel), is considerably re- 
duced in the presence of d electrons. Moreover, these 
small xc effects partially compensate each other, leading 
to an overall effect of no more than 3%. 

Now we focus on the role that localized d bands play in 
the decay mechanism of hot electrons in Cu. This issue 
has been investigated before^ by replacing the various 
contributions to Eq. i|15f) by the corresponding FEG con- 
tributions. Here, we follow a more meaningful procedure. 
As in Rcf . we use Eq. I|15|l (hence, with inclusion of nei- 
ther XC effects nor crystalline local- field corrections) , but 
now we replace the various contributions to the full cal- 
culation by the corresponding s-contributions where the 
3d shell is assigned to the core. These contributions are: 
(i) the imaginary part of the dielectric matrix eG,G(q, ^), 
which represents a measure of both the number of states 
available for the creation of e-h pairs and the coupling 
between states below and above the Fermi level, (ii) the 
denominator |£g,g( c 1j ^)| , which accounts for screening 
effects, and (iii) the matrix elements -Bk.i:k-q./(G) ac- 
counting for the coupling of the hot electron with avail- 
able states above the Fermi level. 

Fig. |2| shows first-principles G°W° calculations of the 
scaled lifetime r(e) x (e — ep ) 2 of hot electrons in Cu, as 
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obtained from Eq. I|15(l . We start with the s-calculation 
(open squares), which we obtain by assigning the 3d 
shell to the core in the generation of the pseudopoten- 
tial. This calculation nearly coincides with the G°W° 
lifetime broadening of hot electrons in a free-gas of 4s 1 
electrons (dashed line), which must be a consequence of 
the fact that band-structure effects are almost entirely 
due to the presence of occupied d bands. In Cu, local- 
ized d bands contribute to the decay of sp hot electrons 
by either opening a <i-band scattering channel at ~ 2 eV 
or by screening the e-e interactions. While the opening 
of the d-band scattering channel modifies the imaginary 
part of the dielectric matrix cg,g(<1, w), the screening 
of d electrons modifies the denominator £g,g(cIj w )| 2 - 
Hence, we have gone beyond the s-calculation (open 
squares) by first including all d states only in the eval- 
uation of Im [cg,g(c[) w)] (triangles) and then including 
all d states also in the evaluation of |eG,G(q, W )P (di- 
amonds). These calculations clearly indicate that: (i) 
occupied d bands yield an slightly enhanced hot-electron 
decay (reduced lifetime) at the opening of the <i-band 
scattering channel (e — Ep > 2eV), which is considerably 
smaller than expected from the greatly enhanced den- 
sity of states at these energies due to the small coupling 
between d states below and sp states above the Fermi 
level for the creation of e-h pairs, and (ii) the key role 
that d electrons play in the hot-electron decay is mainly 
due to screening effects. We also note that the full di- 
agonal G°W° calculation (open circles), which includes 
neither crystalline local-field corrections nor xc effects, 
is very close to the full GWT calculation (solid circles). 
This shows that the overall impact of these effects in this 
material is never larger than a few per cent. 

The screening of d electrons in the noble metals was in- 
cluded by Quinn^ by embedding the ns 1 free electrons 
in a polarizable background of d electrons characterized 
by a dielectric constant instead of unity. The corrected 
lifetime is then found to be larger than in the absence of 

1 /2 

d electrons by roughly a factor of e d . For Cu and Au 
e<2 = 5.6, and for Ag a = 3.4. 3 — Hence, this simple model 
yields corrected lifetimes in Cu that are larger than in a 
FEG by roughly a factor of ~ 2.5, in qualitative agree- 
ment with first-principles calculations. Nevertheless, this 
model cannot account for the existing differences between 
the impact of occupied d states in the lifetime broadening 
of hot electrons in Cu and Au, which both have approx- 
imately the same value of td- 

A comparison between G°W° hot-electron lifetimes in 
Cu and those determined from TR-TPPE experiments 
was presented in Ref. At low hot-electron energies, 
the calculated lifetimes were found to be in agreement 
with the low-energy measurements of Knoesel et al.& At 
larger energies, very good agreement was found with the 
lifetimes measured by Ogawa et al& at the (110) sur- 
face of Cu, the only low-index surface with no band gap 
for electrons emitted in the direction perperdincular to 
the surface. Our results indicate that the inclusion of 
exchange and correlation does not substantially change 
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FIG. 4: As in Fig.Q but for Ag. Here, solid and dashed lines 
represent GWT and G°W° lifetimes of hot electrons in a FEG 
with r s = 3.02, corresponding to the average density of 5s 
electrons in Ag. Inverted triangles represent the TR-TPPE 
measurements reported in Ref. 34. 



this agreement. 



B. Silver 

We have plotted in Fig.0]our full first-principles GWT 
calculation of the average lifetime r(e) of hot electrons 
in Ag (solid circles), as obtained in the ALDA from 
Eqs. lfTUI) - ((T2l but with the inverse dielectric matrix 
Cq G '(l' cj ) re pl ace d by that of Eq. (|19|) . For comparison, 
we also plot in the same figure our full first-principles 
G°W° calculation (open circles), first-principles GWT 
and G a W° s-calculations (solid and open squares), and 
GWT and G°W° FEG calculations with the electron den- 
sity no equal to that of 5s 1 electrons in Ag (solid and 
dashed lines). Also shown in this figure are the TR-2PPE 
measurements reported by Bauer and Aeschlimann^* for 
the lifetime of hot electrons in a thin Ag film (inverted 
triangles). 

Fig.0]shows that both in the presence (solid and open 
circles) and in the absence (solid and open squares) of d 
states the overall impact of exchange and correlation on 
the lifetime of hot electrons in Ag is very small (~ 2%) 
in the whole hot-electron energy range under study. As 
occurs in the case of Cu, the two separate sources of xc 
effects, which are both found to yield significant effects 
in the absence of d states and considerably smaller effects 
in the presence of d states, almost cancel each other. We 
also note from Fig. 0] that at the lowest hot-electron en- 
ergies there is good agreement between our full G°W° 
and GWT calculations and the experimental data, which 
are both ~ 2.5 larger than in the absence of occupied d 
states. At larger energies, however, our calculations lie 
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FIG. 5: As in Fig. |3] but for Ag. Solid circles now represent 
the full GWT matricial calculation of Fig. [I] but now mul- 
tiplied by (e — £f) 2 ■ The dashed line represents the G°W° 
calculation for hot electrons in a FEG with r a = 3.02. 



FIG. 6: As in Fig. Q but for Au. Here, solid and dashed 
lines represent GWT and G°W° lifetimes of hot electrons in 
a FEG with r s — 3.01, corresponding to the average density 
of 6s 1 electrons in Au. 



slightly higher than the experimental curve. 

At this point, it is interesting to notice that while hot- 
electron lifetimes in a FEG with r s = 2.67 (corresponding 
to the average density of 4s 1 electrons in Cu) are larger 
than in a FEG with r s — 3.01 (corresponding to the av- 
erage density of 5s 1 electrons in Ag), hot electrons are 
found to live slightly longer in Ag than in Cu. Both in Cu 
and Ag first-pinciples s-calculations nearly coincide with 
their FEG counterparts, showing that band-structure ef- 
fects are almost entirely due to the presence of d elec- 
trons. Hence, the impact of d electrons is found to be 
larger in Ag (with the onset of the d band at ~ 4 eV be- 
low the Fermi level ejr) than in Cu (with the onset of the 
d band at ~ 2eV below ep). 

Fig. [SI shows first-principles G°W° calculations of the 
scaled lifetime t(e) x [e — £f) 2 of hot electrons in Ag, as 
obtained from Eq. (|15fl by (i) assigning the Ad shell to the 
core in the generation of the pseudopotential (squares) , 
(ii) including all 4d 10 states only in the evaluation of 
Im[eG,G(q,w)] (triangles), (hi) including all Ad 10 states 
also in the evaluation of |eG,G(q, ^)|~ 2 (diamonds), and 
(iv) fully accounting for the presence of 4d 10 states (open 
circles). First of all, we note that it matters very lit- 
tle whether d electrons are included in the evaluation 
of Im [eg,g(<1j u> )] or n °t- This is obviously due to the 
fact that in the case of Ag and at the hot-electron ener- 
gies under study, d electrons are too far below the Fermi 
level to participate in the creation of e-h pairs. Neverthe- 
less, it does not matter how far they are located below 
the Fermi level for them to participate in the screen- 
ing of e-e interactions. Indeed, Fig. [5] shows that the 
screening of d electrons, which enters in the evaluation 
of |eG.G(q, w)p 2 , is responsible for the hot-electron life- 
times in Ag being ~ 2.5 times larger than in the absence 



of occupied d states. 



C. Gold 

Fig. H3 shows our full first-principles GWT calculation 
of the average lifetime r(e) of hot electrons in Au (solid 
circles), as obtained in the ALDA from Eqs. fTU j) -l(T2* jl but 
with the inverse dielectric matrix e G G ,(q, to) replaced by 
that of Eq. I|19fl . For comparison, we also plot in the same 
figure our full first-principles G°W° calculation (open 
circles), first-principles GWT and G°W° s-calculations 
(solid and open squares), and GWT and G°W° FEG 
calculations with the electron density no equal to that of 
6s 1 electrons in Au (solid and dashed lines). The con- 
verged G°W° calculation represented in Fig. HO by open 
circles replaces the lifetimes reported in Ref. 17, which 
were too large by an overall factor of ~ 1.4. 

As in the case of Cu and Ag, the overall impact of ex- 
change and correlation on the lifetime of hot electrons in 
Au is very small, now ~ 2%. We also note that our pseu- 
dopotential G°W° calculations are in excellent agree- 
ment with the all-electron linearized augmented plane 
wave (LAPW) G°W° calculations reported in Ref. H3, 
although they are ~ 20% larger than the correspond- 
ing all-electron linear- muffin-tin-orbital (LMTO) calcu- 
lations reported in Ref. [ijl As noted in Ref. l2lt this 
discrepancy should be attributed to the atomic-sphere 
approximation used in the LMTO calculations of Ref-ITsl 
The LAPW G°W° lifetimes of Ref. were found to ac- 
curately reproduce the BEES spectra for the two pro- 
totypical Au/Si and Pd/Si systems, although they were 
approximately 40% shorter than the TR-2PPE measure- 
ments reported in Ref. ^3- Our GWT calculations indi- 
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FIG. 7: As in Fig. 01 but for Au. Solid circles now represent 
the full GWT matricial calculation of Fig. El but now mul- 
tiplied by (e — £f)' 2 ■ The dashed line represents the G°W° 
calculation for hot electrons in a FEG with r s = 3.01. 



cate that the inclusion of exchange and correlation does 
not change these conclusions. 

A closer analysis of our G°W° calculations is presented 
in Fig. where various calculations of the scaled life- 
time r(e) x (e — Ef) 2 of hot electrons in Au is exhib- 
ited, as obtained from Eq. l(T5|l by (i) assigning the 5d 
shell to the core in the generation of the pseudopotential 
(squares), (ii) including d states only in the evaluation 
of Im [eG,G(q, ^)] (triangles), (iii) including d states also 
in the evaluation of |£g,g(c1, <^)| 2 (diamonds), and (iv) 
including d states everywhere (open circles). As in the 
case of Cu and Ag, these calculations show the major im- 
portance of the dynamical screening of d electrons, which 
is responsible for the hot-electron lifetimes in Au being 
~ 2.5 times larger than in the absence of occupied d 
states. 

As noted in Ref. 0, the role that occupied d states 
play in the screening of e-e interactions is more impor- 
tant in Au than in Cu, although the static polarizable 
background of d electrons in these materials should be 
expected to be the same. However, our calculations indi- 
cate that 5d bands in Au are more free-electron-like than 
3d bands in Cu, thereby allowing Au 5c? electrons to be 
more effective in the dynamical screening of e-e interac- 
tions. 



IV. SUMMARY AND CONCLUSIONS 

We have carried out extensive first-principles calcu- 
lations of the inelastic lifetime of low-energy electrons 
in the noble metals Cu, Ag, and Au, in the framework 
of the GWT approximation of many-body theory. This 
approximation treats on the same footing xc effects be- 



tween pairs of electrons within the Fermi sea (screening 
electrons) and between the excited hot electron and the 
Fermi sea. Our ALDA calculations indicate that the im- 
pact of exchange and correlation on cither screening or 
the bare interaction of the hot electron with the Fermi 
sea, which is significant in the absence of d states, is con- 
siderably reduced in the presence of d states. Moreover, 
these small xc contributions have opposite signs and it 
is the latter which dominates, leading to GWT lifetimes 
that are larger than their G°W° counterparts by ~ 3% 
in Cu and ~ 2% in Ag and Au. 

We have established the role that occupied d states 
play in the relaxation of hot electrons in the noble met- 
als. We have found that deviations from the hot-electron 
lifetimes in these materials from those of hot electrons 
in the corresponding free gas of valence sp electrons are 
mainly due to the participation of occupied d states in 
the screening of e-e interactions, no matter whether oc- 
cupied d states can participate (as occurs in Cu and Au) 
or not (as occurs in Ag) in the creation of e-h pairs. Au 
5c? electrons lie further away from the nuclei than Cu 3d 
electrons, thereby occupied d states in Au being more 
capable to screen the e-e interactions than in Cu. The 
dynamical screening of d electrons yields lifetimes of hot 
electrons that are larger than in the absence of d states 
by a factor of ~ 2 in the case of Cu, and by a factor of 
~ 2.5 in the case of Ag and Au. 

We have found that our G°W° calculations of the 
hot-electron lifetimes in Au are in excellent agreement 
with the corresponding all-electron LAPW calculations 
reported in Ref. 0, which gives us confidence in the ac- 
curacy of our pseudopotential calculations. Our G a W° 
lifetimes are found to be ~ 20% larger than the cor- 
responding all-electron LMTO calculations reported in 
Ref. [l8j, which should be attributed to the atomic-sphere 
approximation used within the LMTO scheme. We also 
note that our G°W° lifetimes of hot electrons in the noble 
metals are systematically higher than those reported in 
Ref. Il6l especially at the lowest energies. These authors 
obtained the lifetime broadening as the full width at half 
maximum of the so-called spectral function, which they 
calculated in the G°W° approximation. Apart from a 
renormalization factor which accounts for the deviation 
of the hot-electron energy from its noninteracting coun- 
terpart and which increases the lifetime by a factor of 
~ 20%, the lifetimes obtained in this way should agree 
with our G°W° calculations. 

Since our GWT calculations are very close to their 
G°W° counterparts, the inclusion of exchange and corre- 
lation does not substantially change the comparison be- 
tween G°W° calculations and TR-2PPE and BEES mea- 
surements existing in the literature for Cu and Au. In the 
case of Ag, we have compared our new G°W° and GWT 
calculations with the recent experimental TR-2PPE mea- 
surements reported in Ref. |3 At the lowest hot-electron 
energies there is good agreement between our calcula- 
tions and the experimental data. At larger energies, our 
calculations are found to lie slightly higher than the ex- 
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